Different thresholds of bond percolation in scale-free networks 
with identical degree sequence 
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Generally, the threshold of percolation in complex networks depends on the underlying structural 
characterization. However, what topological property plays a predominant role is still unknown, 
despite the speculation of some authors that degree distribution is a key ingredient. The purpose 
of this paper is to show that power-law degree distribution itself is not sufficient to characterize the 
threshold of bond percofation in scafe-free networks. To achieve this goal, we first propose a famify 
of scale-free networks with the same degree sequence and obtain by anafyticaf or numerical means 
severaf topofogicaf features of the networks. Then, by making use of the renormafization group 
technique we determine the threshofd of bond percofation in our networks. We find an existence of 
non-zero thresholds and demonstrate that these threshofds can be quite different, which impfies that 
power-law degree distribution does not suffice to characterize the percofation threshofd in scafe-free 
networks. 

PACS numbers: 05.70.Fh, 89. 75. He, 64.60.Ak, 87.23.Ge 
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I. INTRODUCTION 



As one of the best studied problems in statistical 
physics, percolation [ij is nowadays also the subject of 
intense research in the fieid of complex networks 0. 
In a network, if a fraction of its vertices (nodes, sites) 
or edges (links, bonds) is chosen independently with a 
probability p to be "occupied" , it may undergo a per- 
colation phase transition: when p is above a thresh- 
old value Pc, called percolation threshold, the network 
possesses a giant component consisting of a finite frac- 
tion of interconnected nodes; otherwise, the giant com- 
ponent disappears and all nodes disintegrate into small 
clusters. So far, percolation in complex networks has re- 
ceived considerable attention in the community of statis- 
tical physics Q , because it is not only of high theoretical 
interest, but also relevant to many aspects of networks, 
including network security 0, H, [1) 0] ' disease spread on 
networks [8, ^ 0, HH, [11 , etc. 

Since global physical properties of random media alter 
substantially at the percolation threshold, which is cen- 
tral to understanding and applying this process, thus the 
precise knowledge of percolation threshold is extremely 
important p^ . The issue of determining or calculat- 
ing the percolation threshold has been the subject of 
intense study since the introduction of the model over 
half a century ago [ij, [l5|. Despite decades of effort, 
there is stiff no general method for computing the per- 
colation threshold of arbitrary graphs, and rigorous solu- 
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tion for percolation threshold is confined to some special 
cases iM. 3 [13, [13 , such as the Barabasi- Albert (BA) 
network [l9|, two dimensionaf lattice, and some other lat- 
tices. In most cases (e.g. lattices in three dimensions or 
above), the percolation threshold is estimated with nu- 
merical simulations, which are often time-consuming }20i]. 
Thus, finding the threshold exactly is essential to investi- 
gating the percolation problem on a particular graph [T^ . 

Perhaps the main reason for studying percolation in 
complex networks is to understand how the percolation 
properties are influenced by underlying topological struc- 
ture. It has been established that degree distribution has 
a qualitative impact on the percolation. Recent stud- 
ies indicated that in uncorrelated scale-free networks the 
percolation threshold is absent [Bl, [6| . Then a lot of other 
jobs followed, studying the infiuences of other proper- 
ties on the percolation properties in scale-free networks; 
these include degree correlations [13, HH , clustering co- 
efficient nil, and so forth. It was found that, degree 
correlations and clustering coefficient can strongfy affect 
some percofation properties, but they cannot restore a 
finite percofation threshofd in scafe-free networks. This 
raises the question as to whether scafe-free degree distri- 
bution is the onfy ingredient responsible for the absence 
of the percofation threshofd. In other words, whether 
power-faw degree distribution suffices to characterize the 
zero percofation threshold in scale- free networks. 

In this paper, we study the effects of power-law degree 
distribution on the percolation threshold in scale-free net- 
works. To this end, we first construct a class of scale-free 
networks with identical degree sequence by introducing 
a control parameter q. We then study analytically or 
numerically the topological features of the networks and 
show that this class of networks has unique topologies. 
Finally, using the renormalization-group theory, we in- 



2 



vestigate analytically the bond percolation problem in 
the considered networks, and find the existence of non- 
zero percolation thresholds depending on parameter q. 
Our findings indicate that the degree distribution by it- 
self is not enough to characterize the percolation thresh- 
olds in scale-free networks. On the other hand, since our 
networks have the same degree sequence and thus the 
same degree distribution, the model proposed here can 
serve as a useful tool (substrate model) to check the im- 
pact of power-law degree distribution on the dynamical 
processes taking place on top of scale-free networks. 

II. NETWORK CONSTRUCTION AND 
STRUCTURAL CHARACTERISTICS 

In this section, we study the construction and struc- 
tural properties of the networks under consideration, 
with focus on degree distribution, clustering coefficient, 
and average path length ( APL) . 




FIG. 1: Iterative construction method of the networks. Each 
link is replaced by either of the connected clusters on the 
right-hand side of arrows with a certain probability, where 
black squares represent new vertices. 



A. Construction algorithm 

The proposed networks (graphs) are constructed in an 
iterative way as shown in Fig. [TJ Let Ht {t > 0) denote 
the networks after t iterations. Then the networks are 
built in the following way: for t = 0, the initial network 
Hq is two nodes connected by an edge. For t > 1, Ht 
is obtained from Ht-i- We replace each existing link in 
Ht-i either by a connected cluster of links on the top 
right of Fig. [1] with probability g, or by the connected 
cluster on the bottom right with complementary prob- 
ability 1 — The growing process is repeated t times, 
with the graphs obtained in the limit t — > oo. Figures [2] 
and [3] show the growth process of two networks for two 
limiting cases of g = and g = 1, respectively. 

Now we compute some related quantities such as the 
number of total nodes and edges in Ht, called network 
order and size, respectively. Let Ly{t) be the number of 
nodes generated at step i, and Et the total number of 
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FIG. 2; (Color online) Illustration of the first four evolution 
steps of the network growth process for the particular case 
q = l. 
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FIG. 3; Sketch of the iteration process of the network for the 
particular case of g = 0. 

edges present at step t. Then L.u(O) — 2 and = 1. 
By construction (see Fig. [T]), we have Et = 4_E(_i = 4* 
(t > 0). On the other hand, each existing edge at a given 
step will yield two new nodes at next step, this leads to 
L„(t) = 2Et^i = 2 X 4*^1 (t > 1). Then the number of 
total nodes Nt present at step t is 

7V,= ^L„(i,) = |(4* + 2). (1) 

The average node degree after t iterations is {k)t — — 
|rp2, which approaches 3 for large t. 

B. Degree distribution 

When a new node i is added to the networks at a cer- 
tain step ti [ti > 1), it has a degree of 2. We denote by 
ki{t) the degree of node i at time t. By construction, the 
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degree ki{t) evolves with time as fci(t) = 2ki{t — 1) = 
2t+i-ti_ xhat is to say, the degree of node i increases by 
a factor 2 at each time step. Thus, the degree spectrum 
of the networks is discrete. In Ht all possible degrees 
of nodes is 2, 2^ 2^, . . ., 2*"^, 2*; and the number of 
nodes with degree k ~ 2*+^^™ is Uk — Lv{m) — 4™^^. 
Therefore, all class of networks Ht have the same degree 
sequence (thus the same degree distribution) in the full 
range of q. 

Since the degree spectrum of the networks is not con- 
tinuous, it follows that the cumulative degree distribu- 
tion 5il is given by Pcum(fc) = where Nt,k = 
X]jt'>fe "^fe' is the number of nodes whose degree is not less 
than k. When t is large enough, we find Pcumik) ~ k~^. 
So the degree distribution P{k) of the networks follows 
a power-law form P{k) ^ k~^ with the exponent 7 = 3, 
independent of q. Notice that the same degree exponent 
has been obtained in the famous BA network 19]. 




FIG. 4: Second construction method of the network for q = 
1 case that highlights self-similarity: The graph after t + 1 
construction steps, Ht+i, is composed of four copies of Ht 
denoted as {6 — 1,2,3,4), which are connected to one 
another as above. 



C. Clustering CoefHcient 

By definition, the clustering coefficient [13] of a node 
i with degree ki is given by Ci = 2ei/[ki(ki — 1)], where 
ei is the number of existing triangles attached to node 
j, and ki{ki — l)/2 is the total number of possible trian- 
gles including i. The clustering coefficient of the whole 
network is the average over all individual C[s. By con- 
struction, there are no triangles in Ht, so the clustering 
coefficient of every node and their average value in Ht 
are both zero. 

D. Average path length 

Let dij represent the shortest path length from node i 
to j, then the average path length dt of Ht is defined as 
the mean of dij over all couples of nodes in the network: 

dt = — , , (2) 

Nt{Nt ~ l)/2 ^ ' 

where 

denotes the sum of the shortest path length between 
two nodes over all pairs. 

For general q, it is difficult to derive a closed formula 
for the APL dt of Ht . But for two limiting cases of q 1 
and q = 0, both the networks are deterministic, we can 
obtain the analytic solutions for APL. 

1. Case of q — 1 

In the special case (see Fig. [5]), the networks are re- 
duced to the (l,3)-flower proposed in [2^. This limiting 



case of network has a self-similar structure that allows 
one to calculate dt analytically. The self-similar struc- 
ture is obvious from an equivalent network construction 
method: to obtain Ht+i, one can make four copies of Ht 
and join them in the hub nodes. As shown in Fig. 21 
network Ht+i may be obtained by the juxtaposition of 
four copies of Ht , which are consecutively labeled as H^ , 
Hf, Ht, and Hf. Then we can write the sum Dt+i as 

A+i=4A + At, (4) 

where At is the sum of length over all shortest paths 
whose end points arc not in the same Ht branch. 

The paths that contribute to At must all go through 
at least one of the four connecting nodes (i.e., W, X, Y 
and Z in Fig. [5]) at which the different Ht branches are 
connected. The analytical expression for At, called the 
length of crossing paths, is found below. 

Denote A"'^ as the sum of length for all shortest paths 
with end points in H" and Ht , respectively. If H^ and 
Hf meet at a connecting node, A"'^ rules out the paths 
where either end point is that shared connecting node. 
For example, each path contributed to Aj'^ should not 
end at node W. If H" and H^ do not meet. A"''' ex- 
cludes the paths where either end point is any connecting 
node. For instance, each path contributed to A]'^ should 
not end at nodes W, X, Y or Z. Then the total sum At 
is 

At = At^'^ + A]'^ + A]-^ + A?'^ + A?'^ -I- Af ^ - 4. (5) 

The last term at the end compensates for the overcount- 
ing of certain paths: the shortest jjath between W and Y, 
with length 2, is included in A^' and Aj'^; the shortest 
path between X and Z, also with length 2, is included in 
both At^^ and At^^ 

By symmetry, A]'^ = A]''^ = A^^ = A?'^ and A]'"^ = 
At''', so that 

At = 4A^^ + 2At^'^ - 4. (6) 
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In order to find A^'^ and A^'^, we define quantity St as 



St 



i<£Ht 



■iW ■ 



(7) 



Considering the self-similar network structure, we can 
easily know that at time t -\- 1, the quantity St+i evolves 
recursively as 

st+i = 2st + [st + {Nt-l)] + [st + {Nt~l)-2] 

= 4s, + ^(4*-l). (8) 



Using si =4, we have 
1 



St = - (4 + 5 X 4* + 3t X 4*) 
9 



(9) 



On the other hand, by definition given above, we have 



E (d'tvv + djw) 
(Nt - 1) E '^^^ + - 1) E "^1^ 
2{Nt - 1) E ^'W^ 



i=iW 



= 2{Nt-l)st. 
Continue analogously. 



(10) 



A 



t' — E/ 



E (^i^ ~^ dy/X + djx) 



diw<diz 



E i^iZ + dzY + d-jv) 



diz <diw 
j£Hlj^X,Y 



2 E id^w+dwx+dJx), (11) 



i£Hl,iyiW,Z 

diw<diz 
3eH^,3^X,Y 

where the symmetry property has been used. After sim- 
ple calculations, we obtain 



A]'' = 2{Nt - 1) St + {Nt-l)'-2{^ 



2{st + Nt-3)-l. (12) 



Substituting Eqs. ^ and ^ into Eq. we obtain 
after simplification 



At = - (-2-f 11 X 16* + 6t X 16*) . (13) 
9 



Thus 



Dt+i =4Dt + - (-2 + 11 X 16* + 6i x 16*) . (14) 
9 



Using _Di = 8, Eq. (fT^ is solved inductively, 
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16 X 4* + 3 X 16* + 6t X 16*) 



(15) 



Inserting Eq. (|15p into Eq. 0, one can obtain the ana- 
lytical expression for dt: 



, 2 
dt — - X 
3 



16 X 4* + 3 X 16* + 6t X 16* 
4 X 16* + 10 X 4* +4 



(16) 



which approximates t in the infinite t, implying that the 
APL shows a logarithmic scaling with network order. 
Therefore, in the case of g = 1, the network exhibits a 
small-world behavior. We have checked our analytic re- 
sult against numerical calculations for different network 
order up to t = 10 which corresponds to A^io = 1048577. 
In all the cases we obtain a complete agreement between 
our theoretical formula and the results of numerical in- 
vestigation. 



2. Case of q = 

For this particular case, our networks turn out to be 
the hierarchical lattice introduced in [26], which is also 
self-similar. Using a method similar to but a little dif- 
ferent from that applied in preceding subsection, we can 
compute analytically the average path length dt- We 
omit the calculation process and give only the exact ex- 
pression as below: 

_ 22 X 2* X 16* + 8*(21t + 42) + 27 x 4* + 98 x 2* 

* ~ 42 X 16* + 105 x 4* + 42 ' 

(17) 

Equation ([TT]) recovers the previously obtained result 
in [27j and has been confirmed by extensive numerical 
simulations. In the large t limit, dt ~ ^2*. On the other 
hand, for large t limit, Nt ~ 4*, so dt grows as a square 
root of the number of network nodes. Thus, in the case 
of q = 0, the network exhibits a 'large-world' behavior of 
typical node-node distances. 



3. Case of < q < 1 

For < (7 < 1, in order to obtain the variation of 
the average path length with the parameter 5, we have 
performed extensive numerical simulations for different 
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FIG. 5: Graph of the dependence of the average path length 
on the tunable parameter q. 



q between and 1. Simulations were performed for net- 
work with order 10924, averaging over 20 network 
samples for each value of q. In Fig. [51 we plot the average 
path length as a function of q. We observe that, when 
q increases from to 1, the average path length drops 
drastically from a very high value to a small one, which 
shows that there is a crossover between small-world and 
'large-world'. This behavior is similar to that in the WS 
model [Hi. 



III. THRESHOLD OF BOND PERCOLATION 

As discussed in previous section, the networks exhibit 
many interesting properties, i.e., they have the same de- 
gree sequence independent of parameter q; they are scale- 
free and non-clustered; and they display a crossover be- 
tween "large-world" and small-world. All these features 
are not shared simultaneously by any previously reported 
networks. Hence, it is worthwhile to investigate the pro- 
cesses taking place upon the model to find the different 
impact on dynamic precesses compared with other net- 
works such as the BA network. In what follows we will 
study bond percolation, which is one of the most impor- 
tant issues in statistical physics [l[ . 

In bond percolation every bond (link or edge) on a 
specified graph is independently either "occupied" with 
probability A, or not with the complementary probabil- 
ity 1 — A. In our case the percolation problem can be 
solved usirig the real-space renormalization group tech- 
nique d^, [2§, [s^, [3l|, [111, giving exact solution for the 
interesting quantity of percolation threshold. Let us de- 
scribe the procedure in application to the network consid- 
ered. Assuming that the network growth stops at a time 
step t ^ oo, when the network is spoiled in the following 
way: for a link present in the undamaged network, with 
the probability A we retain it in the damaged network. 
Then we invert the transformation in Fig. 1 and define 



n = t — T for this inverted transformation, which is actu- 
ally a decimation procedure [sj. Further, we introduce 
the probability A„ that if two nodes are connected in the 
undamaged network at t = t — n, then at the nth step 
of the decimation for the damaged network, there exists 
a path between these vertices. Here, Aq — A. We can 
easily obtain the following recursion relation for A„ 

A„+i = q (A„ + Xl - AfJ + (1 - q){2Xl - (18) 



fixed points) , among 
•'^ is invalid, 



Equation (flS)) has four roots (i.e. 
which the root A = — — ^ yj^ — &q + q" 
because it is less than 0. The other three fixed points 
are as follows: two stable fixed points at A = and 
A = 1, and an unstable fixed point at Ac that is the 
percolation threshold. The reason for the unstable fixed 
point corresponding to the threshold is as follows: at any 
< Aq < Ac, A„ approaches as n ^ oo, which means 
there is no percolation; while at any Ac < Aq < 1, A„ 
approach 1, indicating an existence of the percolating 
cluster. 

The exact expression of Ac as a function of q is 

Ac = -^ + ^V5-6<7 + g2. (19) 

Interestingly, for the case of g = 0, Ac is equal to , 

which is the inverse of the golden ratio (/) {(f> = ) 
and is the same value as the si te p ercolation threshold 
for the "B" lattice discussed in [3 [i3| ■ We present the 
dependence of Ac on q in Fig. [51 which indicates that the 
threshold Ac decreases as q increases. When q grows from 
to 1, Ac decreases from « 0.618 to 0. 




FIG. 6; The dependence relation of percolation threshold Ac 
on the parameter q. 

Thus, in a large range of parameter q (i.e., q < 1), 
there exists a critical non-zero percolation threshold Ac 
such that for A > Ac a giant component appears spanning 
the entire network, for A < Ac there are only isolated 
small clusters. The existence of percolation thresholds in 
our networks is in sharp contrast with the null threshold 
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found in a wide r ang e of previously studied scale-free 
networks 0, i, M, ■ 

Note that since the susceptible-infected-removed (SIR) 
model can be mapped to the bond percolation prob- 
lem 1^, for the SIR model on our networks the 
epidemic prevalence undergoes a phase transition at a 
finite threshold Ac of the transmission probability. If in- 
fection rate is above Ac, the disease spreads and infects a 
finite fraction of the population. On the contrary, when 
infection rate is below Ac, the total number of infected in- 
dividuals is infinitesimally small in the limit of very large 
populations. The existence of epidemic thresholds in the 
present networks is compared to the result for some other 
scale-free networks, where arbitrarily small infection rate 
shows finite prevalence [33l |. 

From Eq. one can see that for different q, the 

networks have distinct percolation thresholds. As known 
from preceding section, the whole class of the networks 
exhibits identical degree sequence (power-law degree dis- 
tribution) and (zero) clustering coefficient, which shows 
that degree distribution and clustering coefficient are not 
sufficient to characterize the threshold of bond percola- 
tion in scale-free networks. One may ask why the consid- 
ered networks have disparate percolation thresholds. We 
speculate that the diverse thresholds in our networks lie 
with the average path length, which needs further con- 
firmation in the future. 



IV. CONCLUSIONS 

We have demonstrated that power-law degree distribu- 
tion alone does not suffice to characterize the percolation 
threshold on scale-free networks under bond percolation. 
To this end, by introducing a parameter g, we have pre- 
sented a family of scale-free networks with the same de- 
gree sequence and (zero) clustering coefficient. We pro- 
vided a detailed analysis of the topological features and 
showed that the model exhibits a rich structural behav- 
ior. In particular, using a renormalization method, we 
have derived an exact analytic expression for the thresh- 
olds of bond percolation in our networks. We found that 
finite thresholds are recovered for our networks in the 
case of g < 1, which is in contrast to the conventional 
wisdom that null percolation threshold is an intrinsic na- 
ture of scale-free networks. Therefore, care should be 
needed when making general statements about the per- 



colation problem in scale-free networks. 

It should be mentioned that the model generation of 
scale-free networks with the same degree sequence is a 
very common problem in complex network research. Ac- 
tually, in the study of the impacts of other characteris- 
tics (besides degree distribution) of scale-free networks 
on the dynamical processes defined on the networks, the 
interference of power-law degree distribution should be 
avoided. In this case, such a model is necessitated. Tra- 
ditionally, the interchanging algorithm (through rewiring 
two links between^four end points) is frequently used to 
achieve this goal [34] . But this algorithm may lead to dis- 
connection of the whole network. We have shown that the 
scale-free networks proposed here have identical degree 
sequence and are always connected. So our networks can 
overcome above deficiency. They may be helpful for in- 
vestigating how other features (say, average path length), 
other than power-law degree distribution, are relevant to 
the performance of scale-free networks. 

Finally, we stress that since we were only concerned 
with the percolation phase transition point, we merely 
gave the exact position of the percolation thresholds, 
omitting some other properties of bond percolation, such 
as the value of the critical exponents governing behavior 
close to the transition, the complete distribution of the 
cluster sizes, and closed-form expressions for the mean 
and variance of the distribution. All these are worth 
studying further in the future, which is beyond the scope 
of this paper. 

Note added. — A relevant publication [35| about bond 
percolation has come to our attention, where the authors 
showed that different percolation thresholds exist for dif- 
ferent networks having the same degree distribution (not 
degree sequence as addressed in this present paper). 
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